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We investigate the statistical orthogonality catastrophe (StOC) in single-particle and many-body 
localized systems by studying the response of the many-body ground state to a local quench. Using 
scaling arguments and exact numerical calculations, we establish that the StOC gives rise to a 
wave function overlap between the pre- and post-quench ground states that has an exponential 
decay with the system size, in sharp contrast to the well-known power law Anderson orthogonality 
catastrophe in metallic systems. This exponential decay arises from a statistical charge transfer 
process where a particle can be effectively “transported” to an arbitrary lattice site. In a many- 
body localized phase, this non-local transport and the associated exponential StOC phenomenon 
persist in the presence of interactions. We study experimental consequences of the exponential StOC 
on Loschmidt echo and spectral function, establishing that this phenomenon should be observable 
in cold atomic experiments through Ramsey interference and radio-frequency spectroscopy. 


The response of many-body states to local quenches 
can be remarkably drastic; even an arbitrarily weak local 
perturbation can substantially modify the structure of 
the final ground state. As shown in the seminal work by 
Anderson [ 1 ], the overlap between the ground states of 
a metallic (i.e. extended states) system with (|G")) and 
without (IG)) a local perturbation has a power-law de¬ 
cay in the system size F = |(G|G')| ~ A”’*', where 7 > 0 
is a function of the scattering phase shift created by the 
local quench, and thus the overlap vanishes in the ther¬ 
modynamic limit (and hence an infrared ‘catastrophe’). 
This phenomenon is the celebrated Anderson orthogonal¬ 
ity catastrophe (AOC) [2], which has far-reaching conse¬ 
quences in various contexts, ranging from the x-ray edge 
singularity [2] and Kondo effect [3-7] to Luttinger liquids 
[ 8 - 12 ] and resonant tunneling in mesoscopic structures 
[13-18]. Moreover, a power-law multifractal orthogonal¬ 
ity catastrophe (OC) can even occur at random singlet 
quantum critical points [19]. 

Quantum localized systems, both single-particle [20] 
and many-body localized (MBL) [21-24], show a cornu¬ 
copia of intriguing properties, such as violation of ergod- 
icity and emergent integrability [25-27]. MBL has at¬ 
tracted tremendous attention recently due to the break¬ 
down of conventional quantum statistical mechanics and 
the existence of inhnite-temperature dynamical quantum 
phase transitions [28] that appear to defy any description 
within a conventional formalism [29, 30]. Despite exten¬ 
sive research in the subject, the orthogonality catastro¬ 
phe, a compelling phenomenon well studied in diffusive 
metals [2, 31, 32], remains largely unexplored for local¬ 
ized systems [33]. Very recently, in an important pa¬ 
per [34] , a new idea of statistical orthogonality catastro¬ 
phe (StOC) has been introduced for the one dimensional 
(ID) localized Anderson model. 

Here, we study the StOC in single-particle and many- 
body localized systems (see Fig. 1). Through scaling ar¬ 
guments and exact diagonalization, we show that StOC 


gives rise to a wave function overlap that is exponentially 
suppressed in the system size and is much stronger, al¬ 
beit statistical in nature, than the well-studied power- 
law AOC in extended systems. In particular, by study¬ 
ing both ID Anderson and Aubry-Andre (AA) models 
(both with and without interactions) we show that the 
exponential StOC is a generic many-body phenomenon 
in localized systems (both single-particle and many-body 
localized) with a typical wave function overlap Ftyp = 
exp (log F) going as 

Ftyp -exp(-/3A), (1) 

as shown in Fig. 2. Here, /3 > 0 is related to the lo¬ 
calization length and the strength of the local quench. 
Moreover, by using the kernel polynomial method [35], 
we calculate the spectral function and Loschmidt echo, 
which allows us to demonstrate that they manifest cer¬ 
tain unique features that are related to exponential StOC 
and can be directly observed in a cold atom experiment 
through Ramsey interference or radio-frequency spec¬ 
troscopy [36]. Since the fundamental nature of Ander¬ 
son and AA localizations are qualitatively different, with 
quantum interference and quantum gap spectrum playing 
the key role respectively, we hypothesize that all localized 
systems (single-particle or many-body) are characterized 
by the exponential StOC established in the current work. 

Model system .— In this work, we consider a class of 
ID interacting fermionic initial Hamiltonians dehned as 

L/2 

Hi = ^ — J(a]aj+i -|- h.c.) -\- fijUj llujUj+i, ( 2 ) 

i=-L/2+l 

where Oj (aj) are fermionic annihilation (creation) opera¬ 
tors, Uj = ajoj is the corresponding fermion number op¬ 
erator, J is the nearest-neighbour hopping strength, and 
t/ > 0 is the nearest-neighbor interaction. Depending on 
how we choose the local chemical potential, the model in 
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FIG. 1: Nonlocal charge transfer and StOC for the AA 
model, (a) Scaling of the radius of zone of disturbance {Rd) 
with system size in the non-interacting AA model. A lin¬ 
ear scaling Rd ~ 0.031/ is obtained. To suppress statistical 
errors, we averaged over 10"^ different values of cj) (evenly dis¬ 
tributed in [—TT, tt]). Inset: nonlocal adiabatic charge transfer 
for a fixed ((> = Srr/d. A particle on site 0 is transferred to a 
distant location (near site 60) by the time-dependent local po¬ 
tential Vq. (b) Many-body ground state overlaps sorted over 
10^ different values of (j> for the AA model with and without 
interaction. Depending on whether there is a charge transfer 
or not during the adiabatic process, the overlap F is vanishing 
or close to unity, leading to a ‘statistical’ version of OC with 
probability P around 0.25. Here we hx vo = 0.4A = 1.6J. 


Eq. (2) can be reduced to either the Anderson model [20] 
if the local potential fXi is drawn from a uniform random 
distribution between [—FF, VF], or the A A model [37, 38] 
if /j,j = Acos(27raj -I- i^), where A is the strength of the 
potential, with a an irrational number and (j) a random 
phase. Both Anderson and AA models have been exten¬ 
sively studied in the context of many-body localization 
(see Ref. [28, 39, 40]). To unify the notation between 
the two models we will use g to denote W and A re¬ 
spectively for the ID Anderson and AA model. In both 
cases, we start with a localized many-body ground state 
jC?), we then quench the model by introducing a local im¬ 
purity into the Hamiltonian through an onsite potential 
Vo = voalao (vg > 0) at site 0 and arrive at a final Hamil¬ 
tonian I/jp = H/ + Vo, with a new ground state JG"). We 
investigate the asymptotic behavior of the many-body 
ground state overlap (fidelity) defined by F = KGJG')], 
as a function of the system size. 

Charge transfer — We consider a time-dependent 
Hamiltonian H{t) = Hj -|- Vo{t/T), where Hj is a static 
Hamiltonian describing a localized system and Voit/r) a 
time-dependent local potential at site 0, and r is cho¬ 
sen large enough so that adiabaticity is maintained. The 
time dependence of the local potential is Vo(~oo) = 0 
and Vb(oo) = uoOQao. In order to characterize the charge 
response, we define the change of density as Sp{x) = 
{'^x}t—oo {'^x'}t——oo, where {nx}t — (^(^)l^a:|G(t)), 

\G{t = —oo)) = jG) and \G{t = oo)) = jG') denote 
respectively the instantaneous initial and final ground 
states of H{t) [34]. As a result of the local potential, for 


certain disorder realizations such that go < 0, P-o+vo > 0 
charge will be transferred from site 0 to an arbitrary site 
Rf (see the inset of Fig. 1 (a)), with 1 < jRp’j < L/2 
(see our coordinate choice in Eq. (2)), and therefore the 
initial and final ground state wave functions will be or¬ 
thogonal (F —?/ 0 in the thermodynamic limit). Whereas 
for other realizations no charge transfer will take place 
and therefore F « 1. It is in this sense that the StOC 
in localized systems is statistical, depending on specific 
random potential realizations. This is captured in our 
numerical calculations for the AA model with and with¬ 
out interactions as shown in Fig. 1(b), where we plot F 
sorted in magnitude over 10 ^ different random realiza¬ 
tions of (f). For {7 = 0, we find an OC with probability 
P{—Vo < ^0 < 0) ~ 0.25, whereas for U/J=l the prob¬ 
ability is reduced slightly. We also establish StOC for the 
interacting Anderson model [41]. 

In order to quantify the distance the charge will move, 
we consider the so-called radius of the “zone of distur¬ 
bance” over which charge transfer takes place 

\x6p{x)\dx 

Rd = (3) 

Ll /2 \Sp{x)\dx 

As a result of the charge transfer landing on any arbi¬ 
trary site the disorder average oi Rd can have a very 
wide spatial distribution, giving rise to the scaling Rd ^ 
r]{g, U)L, as shown in Fig. 1 (a) for the non-interacting 
AA model, in agreement with Ref. [34] for the ID non¬ 
interacting Anderson model. In addition, we also estab¬ 
lish Rd ~ F in the MBL phase of both interacting ID 
Anderson and AA models [41] . Thus, we have established 
Rd ~ L for the ground state of generic localized systems, 
irrespective of single-particle or many-body localized. 

Finite size scaling of F — Here, we establish the lead¬ 
ing finite size scaling of F in the non-interacting limit 
based on the dynamic charge transfer picture presented 
above. We then generalize the scenario to the MBL 
case. We begin by considering N particles in the lo¬ 
calized phase. The many-body ground state is con¬ 
structed by filling the N lowest single-particle eigen¬ 
states, which have localized wave functions going as 
4>n{x) OC exp(— jcc—x„j/.^), with Xn the localization center 
and ^ the single-particle localization length. Focusing on 
a disorder realization that transfers charge from site 0 to 
site Rf, the leading contribution to the overlap goes as 

F ~ y dx(j)o{x)4>'fix) ~ (4) 

with (f'fiL referring to the post-quench single-particle 
wave function. Consequently, the disorder averaged 
logF ^ —\Rf\/^, where \Rf\ is the averaged distance 
of charge transfer, which has the same scaling as Rd, 
i.e., \Rf\ ~ Rd ~ L. Thus, we arrive at one of our main 
results log F ~ —/3F, and 

Ftyp - exp (-/3F), P = F/4^, 


(5) 
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FIG. 2: Exponential StOC in single-particle and many-body localized systems, (a) Anderson model with vg — OAW = 4J 
and ?7 = 0. log(F’) is obtained by averaging over 10® disorder realizations. The typical fidelity has an exponential decay 
Ftyp ~ i,-0 0724L^ Inset: AOC for a clean (metallic) system with VF = 0, no = 0.2J and [7 = 0. The overlap has a power-law 
decay F ~ t^-sxio ^ expected from Anderson [1]. (b) AA model with vg = 0.4A = 1.6J and U — 0. log{F) is obtained by 
averaging over 10® different values of 4>, which are uniformly distributed in [—rr, tt]. The typical fidelity has an exponential decay 
ftyp ~ e . Inset: OC with power-decaying overlap Ftyp ~ 7 ;, 2 xio delocalized region with A = 0.25no = O.IJ. (c) 

Anderson model a MBL phase with U = J and vg = 0.4IF. log(f’) is obtained by averaging over 10^ disorder realizations. The 
typical fidelity has a exponential decay Ftyp ~ g-o osisi ^ g-o.oi8ii^ £gp ^r _ ^gj 


where P is the probability of creating the charge transfer 
event, i.e. P = vgl{2W) and P{—vg < < 0) for the 

non-interacting Anderson and AA models, respectively. 

The above argument can also carry over to the MBL 
case via the ‘1-bits’ formalism where the ‘1-bits’ are re¬ 
lated to charge degrees of freedom by quasi-local unitary 
transformations and their interactions are exponentially 
local [26, 27]. We then expect from the charge transfer 
picture that the overlap between interacting many-body 
ground states before and after the charge transfer to go 
as F ~ exp(-|i?/|/^MBL), where ^mbl is the many-body 
localization length that can be computed through local 
integrals of motion [42]. Consequently for this many- 
body case, Ftyp still approaches zero exponentially as the 
non-interacting case, yet with /3 = Puyto /We now 
turn to exact numerical calculations to put the scaling 
analysis on a solid footing. 

Numerical results .—We now present our numerical re¬ 
sults in order to establish the exponential StOC for lo¬ 
calized systems. We use exact diagonalization to diag¬ 
onalize Hj and Hp to obtain the energy spectrum and 
eigenstates. For the noninteracting case, the many-body 
ground states are constructed by simply filling the single¬ 
particle eigenstates and the overlaps are calculated via 
determinants (see Supplemental Material [41]). For the 
interacting case, the overlaps can be computed directly 
by inner products of |G) and |G') after an exact diagonal¬ 
ization of the many-body hamiltonian. We then calculate 
the wave function overlap (F) averaged over different dis¬ 
order realizations for both Anderson and AA models. 

Our numerical results for the StOC in the Anderson 
and AA models are shown in Fig. 2, which clearly con¬ 
firms our scaling analysis, and establishes the exponen- 



FIG. 3: /3 dependences for the Anderson ((a), (b), and (c)) 
and A A ((d), (e), and (f)) models in the noninteracting limit. 
The other parameters are chosen as vg = 0.4IF (vg = 0.4A) 
in (a) and (b) ((d) and (e)); W = 10J in (c) and A = 4J in 
(f). Dashed magenta lines are hts to the linear relationships 
in each case. 


tial StOC Ftyp ~ exp(—/3F). For the sake of comparison, 
we have also plotted the results for the clean (metallic, 
FF = 0) limit of the Anderson model and for the delocal¬ 
ized phase of the AA model (A/J < 2), respectively. For 
these extended states, the overlaps follow a power-law de¬ 
cay, as expected from AOC [1]. For the non-interacting 
case we are able to reach sufficiently large system sizes, 
whereas in the MBL phase we are limited to F = 24 as 
our largest system. 

In Fig. 3, we show the dependence of /3 on IF, A,^ 
and vq in the non-interacting limit. First, we focus on 
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the strength of the on site potential, which is related to 
the localization length as 1/^ ~ \g ~ 9cY near with 
a localization length exponent v (for the ID Anderson 
and AA models, gc is 0 and 2A, respectively). Deep 
in the localized phase ^ saturates to the lattice spacing, 
whereas it diverges on approaching the localization tran¬ 
sition point. We thus expect /3 to be zero at g = gc 
and to increase as we increase g from gc, which is in 
excellent agreement with our numerical calculations. In 
addition, the scaling analysis also predicts /3 oc 1/^. We 
compute the localization length directly from the decay 
of the single-particle wavefunction 4>{x) from its maximal 
value (in absolute value) and then average over disorder 
realizations. Close to the localization transition and deep 
within the localized phase such a method of calculating 
^ suffers from large numerical error, nonetheless in the 
intermediate parameter regime we establish a clear lin¬ 
ear relation. Finally, we come to the dependence of /3 on 
vq for a fixed g. The results are in good agreement with 
P (X Vq, which allows us to conclude that our numerical 
results are consistent with /? = P/4^, and thus provide 
another independent confirmation of our scaling analysis. 

Experimental detection. —The AOC has far-reaching 
consequences for several dynamical phenomena, such as 
X-ray edge singularity in metals [2] and power-law sup¬ 
pressed Loschmidt echo for ultracold atoms in the long 
time limit [43]. We expect the exponential StOC in lo¬ 
calized systems to have significant implications in similar 
experimental settings. In a setup proposed to probe the 
OC physics [43-45], a nature experimentally measurable 
quantity to characterize the exponential StOC is the typ¬ 
ical Loschmidt echo defined as 5'typ(t) = exp(log [^'(t)]), 
where 

S{t) = {G\e^^^^e-^^^^\G). (6) 

In Fig. 4 (a) we show Styp{t) for the Anderson model 
[46]. As a result of the StOC and the localized nature of 
the many body ground state we find S'typ(t) saturates to 
a nonzero value in the long-time limit. This is in sharp 
contrast to the metallic case where lS'(t)l has a power- 
law decay due to the AOC [43]. Another quantity that 
is also directly measurable in experiments is the spectral 
function (which is related to the “hole” Green function in 
the AOC literature [2]) A{uj) = J^n' I {G\n')pS{uj—{En' — 
Eq)), where \n') is the n'-th eigenstate of EIp with eigen 
energy En' and Eq is the ground state energy of Etj. 
Using the kernel polynomial method [35], we calculate 
the disorder-averaged spectral function for both the AA 
and Anderson models (see [41]). As shown in Fig. 4(b) 
and (c), the spectral function exhibits two peaks with 
exponentially decaying tails, which are unique features 
of the exponential StOC in localized systems (whereas 
in the AOC, A{uj) has only one peak and a power-law 
dependence [2]). 

With the recent experimental realizations of single¬ 
particle [47-49] and many-body localization [50-52], the 



FIG. 4: Numerical simulations of the typical Loschmidt echo 
and spectral functions, (a) The typical Loschmidt echo as a 
function of time for the Anderson model with and without in¬ 
teraction. Inset: Loschmidt echo in the clean limit, which has 
a power-law decay at long time, (b) The disorder-averaged 
spectral function for the interacting Anderson model. Unlike 
the metallic or the ordinary band-insulator cases, A{uj) for 
localized systems has two divergent peaks at either cu = 0 
or uj = Vo and non-zero weight between (0,Uo)- (c) Zooming 
in on the red rectangular region of (b). Near w = 0, an ex¬ 
ponential decay A(uj) ~ 6 x -I- 0.1 is obtained. The 

parameters are chosen the same as in Fig. 2 (see [41] for more 
details). 

above intriguing features could be experimentally ob¬ 
served with present techniques. Considering a ultracold 
atom experiment for instance, the localized ground states 
could be achieved by adiabatic ramping and the local 
quench can be switched on and off by manipulating lo¬ 
cally two-level atoms (the impurity atoms of a strongly 
imbalanced mixtures [53, 54]), and the Loschmidt echo 
and the spectral function can be directly measured 
through Ramsey interference and radio-frequency spec¬ 
troscopy, respectively (see Ref. [43-45] for details). In 
electronic materials, nontrivial effects of the exponential 
StOC may have already been observed in recent exper¬ 
iments of quantum electron glasses [55, 56], where dra¬ 
matic sluggish dynamics have been reported. 

To conclude, we have established a statistical expo¬ 
nential orthogonality catastrophe in both single-particle 
and many-body localized systems. We introduced scaling 
arguments to explain this ‘stronger’ catastrophe. Using 
exact numerical calculations we have verified the scaling 
analysis and its generalization to the MBL case. The ex¬ 
ponential StOC can be used as a diagnostic tool to iden¬ 
tify localization and could be directly observed in cold- 
atomic or other mesoscopic systems with current tech¬ 
nologies. 
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I. SUPPLEMENTAL MATERIAL EOR “EXPONENTIAL ORTHOGONALITY CATASTROPHE IN 
SINGLE-PARTICLE AND MANY-BODY LOCALIZED SYSTEMS” 

A. Calculations of ground-state overlaps for free fermions 

In this section, we provide the details on how we calculate the ground-state overlap for free fermions. As in the 
main text, we denote by |G)and |G') the ground states of Hj and Hp- Since we only consider the free fermion case. 
Hi {Hp) can be written as Hj = J2ij {Hp = J2ij '^liHp)ijaj), where Hj {Tip) is the so called kernel 

Hamiltonian. One can diagonalize Hi {Hp) by a unitary transformation U (V): Hi = U^SiU {Hp = V^EpV) to find 
the single-particle eigenmodes bj = HjkO-k {cj = J2k '^jkO-k)- Here £i = diag(e{, £21 ’ ’ ’) {^f = diag(ef, e^, • • •)) is 
a diagonal matrix. Suppose the lattice size is L and the total particle number is N. For free fermionic systems, the 
ground states of Hi and Hp are constructed by filling the lowest N single-particle eigenstates: 

N 

iG) = 

3 

|G') = 

3 

Thus the overlap reads: 

N N 

= |(G|G')| = |(0|P[&.n4lO)l- (7) 

i k 

Computing F directly from Eq. (7) is challenging. However, we can simplify Eq. (7) by expressing the eigenmodes 
c^s as a combination of b\s c\ = Y.m(y’^)kmaln = which corresponds 

to a change of basis. Defining another matrix A = V^U, we have c\ = ^kub^ and consequently 

N N 

F = i(oin^.nE^'="^nio)i- 

j k n 

Noting that {bj,bl^} = 6jn, bn\0) = 0 and = 0 for all n and p > 1, we arrive at a greatly simplified formula 

F = \detB\, (8) 

where = Aij (1 < i,j < N) is a matrix with elements taken from A. The above equation is very useful in 
our numerical calculations. With this formula, we are able to calculate the overlap efficiently for different random 
realizations. 


B. More results on charge transfer and orthogonality catastrophe 

In the main text, we have established an exponential statistical orthogonality catastrophe in localized systems and 
its relation to adiabatic charge transfer. Here we show additional numerical results to support our claim. 


1. The Anderson model 

In Fig. 5, we plot the scaling of Rp with system size and the StOC for the Anderson model. From Fig. 5(a), 
one indeed obtain a linear scaling Rp ~ L for the interacting Anderson model in the MBL phase, thus validating 
the heuristic arguments in Ref. [34]. From Fig. 5(b), we have statistical orthogonality catastrophe due to the 
nonlocal charge transfer, as expected. The probability of StOC for the noninteracting Anderson model is about 
P « vo/{2W) = 0.2, whereas this probability is reduced for the interacting case, similar to case of the AA model (see 
Fig. I in the main text). 
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FIG. 5: Scaling of Rd and many-body ground state overlaps for the Anderson model, (a) Scaling of Rd with system size 
for the interacting Anderson model. Linear scalings Rd ~ 0.321/ and Rd ~ 0.27L are obtained for IF = 5J and W = lOJ, 
respectively. Here we have fixed U = J and used 10^ random realizations to suppress statistical errors, (b) Many-body ground 
state overlaps sorted over 10^ disorder realizations for the Anderson model. For the interacting (noninteracting) case, the 
lattice size is L = 24 (L = 200). Other parameters are chosen as vo = 0.4tF = 4J. 



FIG. 6: Nonlocal charge transfer in the interacting Aubry-Andre model with U = J. Rd is obtained by averaging over 10"^ 
different values of (j) evenly distributed in [—tt, tt] and vq is chosen as Vg = 0.4A. 


2. The Aubry-Andre model 

In the main text, we have shown the StOC in this model for both interacting and noninteracting cases (see Fig. 
1(c)). Here we show more numerical results on the charge transfer and StOC for the interacting Aubry-Andre model. 
In Fig. 6, we plot the scaling of the radius of zone of disturbance Rjj [Rd is defined by Eq. (4) in the main text) with 
the system size. One may notice that the numerical lines are not exactly straight lines and there are large deviations 
due to the finite-size effect. Yet, it is clear that the qualitative trend is still correct. In Fig. 7, we plot the scaling of 
—log(E) with the system size for the interacting AA model. Here the finite-size effect is also non-negligible. In order 
to clarify the role of finite size effects we compare the many body interacting case with the non-interacting results 
that can reach much larger systems sizes as shown in Fig. 8(a). From this figure, it is clear that the influence of 
interaction on —log(F) is weak deep in the localized region. Despite the step like features at small L, as shown in Fig. 
8 (b) and (c) this is truly just an artifact of small system sizes, where increasing L gives rise to a clear —log(F) ~ L 






















FIG. 7: Exponential statistical OC in the interacting Aubry-Andre model. Due to finite size effect, the numerical data has 
a large deviation from the corresponding fitted line. —log(E) is obtained by averaging over 10^ different valnes of (f) evenly 
distributed in [—tt, tt]. Other parameters are chosen as U = J and vg = 0.4A. 



FIG. 8: (a) Gomparison of exponential statistical OG in the interacting and noninteracting Aubry-Andre model. In the deep 
localized region (A = 10J), the interaction does not modify —log(f’) very much, (b) and (c) plot the scaling of —log(F') with 
system size for the noninteracting case. Results are obtained by averaging over different values of 0 evenly distributed in 
[—TT, tt] and here vg = 0.4A = 4J. 


system size dependence. Thus, as shown in these two figures, the linearity of —log(T’) with L becomes more and more 
evident as the system size becomes larger and larger. As we have shown, deep in the localized phase, the presence 
of interactions does not quantitatively change —log(F) by a large amount, we thus conclude that —log(J^) also goes 
linearly with L, giving the exponential scaling Ftyp ^ as expected. 


C. Calculations of the Loschmidt echo and spectral function 

1. Calculations of the Loschmidt echo 


In the main text, we have calculated the Loschmidt echo defined as: 

S{t) = 


(9) 
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FIG. 9: Disorder-averaged spectral functions, (a) The Anderson model without interaction. The parameters are chosen as 
W — lOJ and vo = 4J, and the spectral function is averaged over 10"* random realizations, (b) The AA model without 
interaction, (c) The AA model with interaction {U = J). Other parameters are chosen as A = 10J and vo = 4J, and A(co) is 
averaged over different values of (j) evenly distributed in [—tt, tt]. For all these three plots, the system size L = 20. 


A direct calculation of S{t) using Eq. 9 confronts lots of computational complexity. For the noninteracting cases, 
the calculation can be greatly simplihed. Note that |G) = ^]|0) is the ground state of Hi, then we have S{t) = 

To calculate (G|e“*^^*|G), one can do a time-evolution of each eigenmode of Hi 




EE 

m k 


( 10 ) 


After bj(t) is obtained, we one can re-express hj{t) as a combination of hi's and then the calculation of S(t) is reduced 
to the computation of a determinant at each time t, following the same procedures in Sec. I. In Fig. 4(a) of the main 
text, we have shown the typical Loschmidt echo Sty-p for the noninteracting Anderson model (the yellow curve). In 
this case, our lattice size is L = 200 and we used 10^ random realizations. For the interacting case, calculations of S{t) 
requires full diagonalizations of many-body Hamiltonians and our system size are limited to L = 14. The numerical 
results of S{t) for the interacting Anderson model are also shown in Fig. 4(a) (the black curve) in the main text. 
For both interacting and noninteracting cases, S(t) has various distinctive features that are related to the exponential 
StOC and are different from that of the Anderson OC for a metallic phase. We have also calculated the Loschmidt 
echo of the AA model for both interacting and noninteracting cases, and very similar features are obtained. 


2. Calculations of the spectral function 


In the main text, we have defined the dynamic spectral function as 

A(a.) = ^ |(G|n0p5(a; - (F„, - Eg)). (11) 

n' 

A direct calculation of A{uj) using Eq. 11 would also require a full diagonalizations of the many-body Hamiltonians 
to obtain the eigenstates \n'). We can avoid this by using the so called kernel polynomial method (KPM) [35] to 
calculate A[uj) for system size up to L = 24. For each random realization of Hj, we use the Lanczos algorithm to find 
the ground state |G) and the corresponding eigen-energy Eq- We then define a new Hamiltonian H' = Hp — Eg to 
shift the eigen energies oi Hp and follow the standard procedures of KPM (see Ref. [35] for details). In order to fit 
the spectrum into the interval [—1,1] (within which the Chebyshev polynomials are defined) we apply a simple linear 
transformation to H'-. 


H' = {H' - b)/a, 


( 12 ) 


where a = — EH^)/{2 — e) and b = (if^ax + ^min)/2i Here e is a small cutoff parameter introduced to avoid 

stability problems (e = 10“^ is chosen in our numerical calculations); EA^ are the maximal and minimal 

eigen energies of H', respectively. Both F(,jax EH^ can be obtained by directly diagonalizing H' using Lanczos 
algorithm. Then the spectral function can be expanded in Chebyshev polynomials [35] 


1 / 

goi^o + 2 E 9kJ^kTk{Cj) 

rv 1 - V 


A{Cj) = 


(13) 
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where rfc(a;) = cos[fc arccos(a;)] are the Chebyshev polynomials and = [{K—k+l) cos{K-\-l) 
are the coefhcients of the Jackson kernel; the moments Vk = {G\Tk{—H')\G) can be obtained iteratively based on the 
recursion relations of the Chebyshev polynomials 

To{x) = 1, T'_i(x) = Ti{x) = X, and Tk+i{x) = 2xTk{x) - Tk-i{x). (14) 

After A{u}) is obtained, one can do a rescaling Cj ^ uj = aCu + 5 to obtain A(ui). In Fig. 4 (b) in the main text, we 
have plotted the disorder-averaged spectral function for the interacting Anderson model. In this plot, we have L = 20 
and used 10^ random realizations. The spectral function is normalized according to f A{uj)duj = 1. As discussed in 
the main text, we have found unique features related to the exponential StOC in localized systems. In Fig. 9, we 
also plot the spectral functions for the noninteracting Anderson and interacting AA models. Distinctive features of 
double peaks and exponentially decaying tails are also obtai ned in these cases. _ 



